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We present a canonical formulation of gravity theories whose Lagrangian is an arbitrary function 
of the Riemann tensor. Our approach allows a unified treatment of various subcases and an easy 
identification of the degrees of freedom of the theory. 



I. INTRODUCTION 

0\ \ . . n • 

, Since H. Weyl introduced them in 1919 pj], theories of gravity whose action is nonlinear in the Riemann tensor 
(contrarily to Hilbert's) have been part of the "landscape" of General Relativity and its various extensions, up to 
_ present-day string theories. In this paper we shall be interested in the action: 

< 



X 



^ JM 

Sg is a functional of the metric g^v{x'') , D is the dimension of the spacetime M. , and / is an arbitrary function of 
the Riemann tensor TZf^^pa ■ ^ 

The Euler-Lagrange equations of motion derived from metric variation of (jl.ip are 



X: 7^(^,.-^-2VpV.— ^-i/g-'' = ^^^ (1.2) 

5_L| Ol<-u)Xpa Ol<-p{p.v)a ^ 

, where is the covariant derivative associated with g^^, , T'^" is the energy-momentum tensor of matter, and f{p^) = 
{fpf + /i^/^)/2 . These equations are generically fourth-order in the derivatives of the metric. To convert them into a 
I [ set of first-order differential equations, one must introduce as new variables some well-chosen functions of the metric 
J> ■ and its derivatives. The identification of these extra degrees of freedom (besides those of Einstein's gravity) is a 
' prerequisite to introducing non- minimal coupling to matter. It is also an important step to study for example the 
. well-posedness of the initial value problem and the number of independent Cauchy data, the global charges associated 
^^O with the solutions, the stability of the theory and the positivity of energy as well as the junction conditions. 

, At the linear approximation around Z? = 4 flat spacetime and in the case f = a Ti? + [3 TZfj^i, TZ^" the identification 
00 ' of the extra dynamical degrees of freedom was made long ago, see 0: They are generally six, corresponding to 
[ a "massive spin-0" together with a "massive spin-2" field (a "ghost" of negative energy), and reduce to five if / 
■ is the square of the Weyl tensor {(3 = —3a). The action was then recast in first-order form and the constraints 
[ were analysed in 0] and [3| with the result, among others, that when / is the square of the Weyl tensor there is 
J> . an additional, "first-class," constraint which generates conformal transformations (see also [1, Q)- The analysis of 
/ = aTZ"^ + (3 TZ^„ TZ^'^ -\- 7 TZ^^^pa TZ^^^p"^ in D dimensions was performed in with the result that the number of extra 
degrees of freedom is reduced if / is the square of the Weyl tensor (a = 2 , /3 = —4 (Z) — 1) , 7 = (_D — 1) (Z? — 2)) or 
' the Gauss-Bonnet combination {a — \ , (i = —4 ,7 = 1). Finally, a canonical formulation of f{TZ) theories has been 
proposed in Q (see also [§|) and the Lovelock case [l3| was treated in [111] . 

Our aim here will be to unify these results and generalise them to the full action (jl.ip . To do so we shall follow 
the procedure of Arnowitt, Deser and Misner (ADM) [l^l and present a canonical formulation of theories of gravity 
yielding the field equations (|1.2p . In contradistinction with previous approaches which consist in choosing as the extra 
variables either the extrinsic curvature of the ADM foliation Q or its time derivative , we shall ascribe a leading 
role to the components TZi^ oc Vtij of the Riemann tensor (in an ADM coordinate system adapted to the foliation). 
This will allow us to write the canonical equations of motion in a compact form. The extra degrees of freedom will 
be encoded in the D {D — l)/2 components of some spatial tensor ^I'*^ and their number will depend on the number 
of components of flij that can be extracted from the equation 

2*'' + J7- = 0- (1-3) 



Our conventions are: TZ^upa = (1/2) {dfipgi/a — di^agp,p) + ■ ■ ■ ; g is the determinant of gi_w ; the signature is ( h + ■ ■ ■ ) ; spacetime 

indices (fi, u, ■ ■ ■) run from to D ~ 1 ; space indices • • ■ ) will run from 1 to D — 1 . 
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The paper is organised as follows. In Section [IT] we introduce a number of auxiliary fields in the action, perform its 
ADM decomposition and simplify it by using some of the constraints which simply determine some of the auxiliary 
fields algebraically. This leads us to our second-order form Lagrangian, see (|2.18p . In Section IIIII we first obtain 
the Hamiltonian of the theory, see (|3.4p and (|3.5p . and, second, write Hamilton's equations of motion, see (|3.9p 
and (jS.lOp . Section IIVI presents the way to reduce the Hamiltonian using constraints and shows how a number of 
well-known subcases are recovered (General Relativity, f(TZ) and "Weyl^" theories). Finally we make the link in an 
Appendix between our formalism and the commonly used Ostrogradsky one. 

We shall present in this paper only the thread of the argument. The details of the calculations will be presented 
elsewhere Il3ll. 



II. 



CHOOSING THE ACTION 



Introducing auxiliary fields 



In order to turn the equations of motion (|1.2p into a set of first-order differential equations we shall start our 
analysis, not with (|l.ip but with the related action 



pcrl 



M 



(2.1) 



The two auxiliary fields g^vpa and ipt^'^P'^ have all the symmetries of TZpi^vpa and are chosen to be independent of each 
other and of y^j, . The reason for introducing them is that the second derivatives of the metric appear only linearly 
in (|2.ip and will be eliminated by means of an integration by parts, see below. If matter does not couple to ipf^'^P"' 
and gfj^i^pa , variation of (|2.ip yields a set of field equations which is equivalent to (|1.2p . Indeed: 



6S 



M 



-Sf^" Sg,,, + in^,p„ - gp,pa) Sifi^'"'"' + 



df 



dgau 



fivpcr 



M 



where 

Ignoring the divergence term for the time being, the field equations are 
and, if the matter action does not depend on g^upa and •p^'^P'^ ,2 
together with 

df 



if 



(2.2) 

(2.3) 

(2.4) 
(2.5) 

(2.6) 



Substituting these extra equations into (j2.4p , we recover the original fourth-order equation of motion ([TT 

A remark is in order here. As is well known fl6i|, one can plug the constraint (|2.6p back into the action (|2.ip and 
consider it as a functional of g^u and g^upa only. Indeed variation of this new action still yields the same equations 
of motion, the only difference being that (|2.5p is replaced by 



dg^upa dgafijs 



{TZafljS — gafljs) — . 



(2.7) 



Whether this equation can be inverted to yield (|2.5p or not imposes an analysis of the various subcases before any 
canonical treatment. For this reason we shall refrain from using (|2.6p straightaway and shall stick to the action (|2.ip . 
This will allow us to treat the different subcases in a unified manner. 



^ For an example of non-minimal coupling to matter, see e.g. Il4l. llSH . 
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B. ADM decomposition 

Suppose that Ai can be foliated by a family of spacelike surfaces St, defined hy t — .'^ Let /ly = gij\xO=t 
with i,j running from 1 to I? — 1 be the metric on , /i its determinant, h^^ its inverse, and Di be the associated 
covariant derivative. Introduce the future-pointing unit normal vector n to the surface , that is, to the basis vector 
fields di with components 5f ; the components of n are rn = , no — —^/y/~g^, "nP — \/—g^, — — ff"Vv^^- 
Decompose then the timelike basis vector do (with components Sq) on n and the spatial vectors di: 5q — N n^^ + fi"^ 5^ . 
N = l/y— 5°° and = —g^^/g^^ are the "lapse" and "shift," respectively; together with the induced metric hij , 
they constitute the "ADM variables" • In terms of these variables we have 

no = -iV, n,; = 0, n° = ^, = (2.8) 
The components of the spacetime metric read 

500 = -N'^ + Pi , goi = A , gij = hi^ , 

and ^—(7 = nVK Introduce finally the extrinsic curvature of Et : 

K,, = V,n, = ^ {k, - i?,/3, - Z?,/?0 , (2.10) 

where a dot denotes a time derivative: /i^j = dthij . 

A standard calculation using (IMl), (US and (see e.g. [13, [11 [11 

for a geometrical derivation) yields the 
Gauss, Codazzi and Ricci equations, that is, the components of the Riemann tensor in terms of the ADM variables: 

'R-ijki = Kik Kji — Kii Kjk + Rijki , 

TZijkn = Tlijkii = DiKjk ~ DjKik , (2.11) 
7^™,n EE 7^,^,, - -N'^ {K,, - £pK,,) + [K ■ K),, + N'^ D,,N , 

fc 



where Rijki is the Riemann tensor of hij , [A ■ B)ij = Aik Bj , Dij = DiDj , and the short-hand notation £0 denotes 
the Lie derivative with respect to the shift: £fiKij = f3'^DkKij + Kik Djl3^ -f Kjk Di(3*' . We can thus perform the 
following decomposition: 



^Mi'P'^ (JZf,i,pa- - Qf_iypcr) = (t>'^^^ (Jltjkl - Pi]kl) - 4 0'^'' {Tlijkn - Pijk) - 2 ^''^ (T^mjn - ^ij) , (2-12) 

where we have introduced the following spatial tensors evaluated on Et : 



Pijkl — Qijkl , Pijk — Qijk^ , ^ij — n^ 71 Qip,ji/ , 

^rjkl ^ ^un ^kp ^Iq ^^^^^ ^ ^rjk ^ ^kn ^^^^^^^ ^ = _2 /l^^' /l^' Tl'^ ^kt^lu 



(2.13) 



In order to pass to canonical formulation, one must remove second derivatives from the action. It is done by 
integrating the second derivative by parts to cast it into a total divergence. A side exercise shows how the time 
derivative of the extrinsic curvature is transformed into a time derivative of '^^^ : 

{k,j - £pK,j) - V^(n^ K-^)~K{K-'^)~N-^ K,, - £p-^'') , (2.14) 

where K = h'^ K^j a.nd A ■ B = Ay B'^ . 

Armed with these preliminaries we can decompose the action (|2.ip as 

5[%,iv,/3^vI/*^^]„,<^^^'=^(/.*^■^p,,fc^,p^,fe] = / d^'xic+dpi^n'^K-^)], (2.15) 

Jm 



^ We here work in an ADM coordinate system. See e.g. [iTI . llSl . Il9| | for a more geometrical approach. 

* Here and in the following indices of (D — l)-dimensional spatial tensors are moved with the induced metric hi 
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with 



(2.16) 



where T^-^fe; and TZijkn are given in (|2.1ip and where it is understood that f{g^iypa) is expressed in terms of pijki , 
Pijk and Qij (see below for examples of such a decomposition). 

As one can see from (|2.16p and (|2.1ip the spatial tensors </)*'''^' and cf)^^'' are not dynamical and their equations of 
motion {dC/ dcf)^^^^ — , dC/dcjf^*' — 0) are simple constraints: 



Pijki — Kik Kji — Ku Kjk + Rijki , Pi]k — DiKjk — D^Kik , 



(2.17) 



which can harmlessly be incorporated into C , see [16l | . As a consequence the Lagrangian density of our theory reduces 
to 



C* = VhN 



(2.18) 



where the spatial tensors ptjki and pijk in f{Qp.vpa) are given by (|2.17p . 

The divergence in (|2.15p is canceled by adding to the action the generalisation of the York-Gibbons-Hawking term 
dSHH (see also ^) : 



S 



dM 



dE^ n'' * • ivT , 



(2.19) 



where dS^ is the normal to the boundary dAi being proportional to the volume element of dAi . Here dS^ is 
outward-pointing if spacelike and inward-pointing if timelike. The field equations (j2.4p . (|2.5p and (|2.6p then derive 
from a Dirichlet variational principle where the induced metric hij and 5'*-' have to be held fixed on dM . 

Let us note that there can be constraints on 5'*-' depending on the form of / (see Section IlII BI for general discussion 
about this constraint) and in some specific cases VE'*-' is given as a function of K^j . Then the surface term (|2.19p cannot 
be used to remove the second derivatives from the action. This is for example what happens in the Lovelock case 
[Tst . For simplicity, we do not consider such exceptional cases in this paper and the divergence in (|2.15p is hereafter 
discarded. More details about the surface term will be presented elsewhere [l3j . 



III. CANONICAL FORMALISM 



A. First-order action and Hamiltonian 



Momenta conjugate to the dynamical variables hij and ^'^-^ are defined as, recalling the definition (|2.10p of Ki. 



n, 



5C* 

6 hi. 



SC* 



Vh 
2 



h"-'-^ -K + K + 2 («- • if )(''^) + N-^ - i"/?**^') - 2 



df 
dpikji 



Kki 



+ N-'Dk N 



df 



dpk 



(3.1) 



= -^/hK, 



where it is understood that the derivatives df /dpijki and df /dpijk have all the symmetries of pijki and pijk , respec- 
tively, and we have used the constraints (|2.17p and discarded a total divergence in p^^ Inversion yields the velocities 



^ When integrated over space this divergence yields the following contribution to the momentum: (1/2) f^dSk {df /dpk{ij)) i where dSk 
is the outward-pointing normal to the (D — 2)-dimensional sphere S being proportional to the volume element of S . Such a term, as well 
as all spatial divergences, will be discarded in this paper but are important when studying, e.g., the energy of the system or junction 
conditions [isll . 
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in terms of the canonical variables (see (|3.9p below for their explicit expression) and the first-order Lagrangian then 



IS 



with 



C* ^p-h + U-ii -H* - d^i^V) 



where 7i* is the Hamiltonian density: 
with 

c = ^ ((* • n) n + • n • n - 2p • n) + 



(3.2) 
(3.3) 
(3.4) 

(3.5) 



where 11 = ft,'-' liij and A • i? • C = y4*j C'^i , and the spatial tensors pijki and py/j in fig^upa) are now given by 

n,fc n,, - n,, n 



life lij; - iii; lljfc / Iljfe \ „ / Ilife 



(3.6) 



B. Hamilton's equations and the algebra of constraints 



The equations of motion consist first in two sets of constraint equations, 67i* /5N — and SH* /5(3'^ ~ , which are 

C = , = . (3.7) 
The equation of motion for Qij , STi* /Sflij = , is also a constraint: 

Of 



an. 



0. 



(3.8) 



Depending on the function fig^i^pa) , this equation may or may not be invertible to give all the components of flij in 
terms of ^''^ as well as hij , 11^ and their spatial derivatives. As announced in the Introduction, the number of extra 
degrees of freedom beyond those of General Relativity will depend on the invertibility properties of p.8p .^ We shall 
come back to this issue at the end of this Section. A systematic way to reduce the action using (|3.8p will be presented 
in Section HVl on some specific examples. 

As for the dynamical equations, the first set, hij ~ STL* jSp^^ and "^^^ = SH* /SHij , gives the velocities in terms of 
the canonical variables: 



w'i = £p^'i + ^ (n f*^' + (* • n) h'^ + 2 (* • n)('^) - 2p'^) - 2 iv 5^ - a (n 

Vh Vh opikji \ clPk{i]) 



(3.9) 



The second set, p'^ = -5H*/5h,j and = -SH* / S'H''' yields, using 



N 



^Vh£p[^] +^ [(* • n) n'^ + ^ki ff* ff' + ft'^ ((* • n) n + * ■ n • n - 2p • n)] 

+ [Dki'f'' d^N - 2 d^^N) + h'^ {NDki-^''^ - -^^^ DkiN)] + — ^^^^ , 
2 2 ohij 

riy ^^/h£p {^j^ ^ ^ + (n • n)y ) + \/ft (iv r!y - d,jN) , 



(3.10) 



Notice that nonlinearity of / in £7ij = 7?.injn is essential for the invertibility of l|3.8| l as well as for the appearance of fourth-order 
time-derivatives in 111.21 1. 
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where the last term of the right-hand side of the first equation reads, using (|3.8 



7 - . ^ ^ khn ^ ~v ^ J-^kl \ ""o / i o — ---^ „ 

Ofty OP'j)klm \ ClPk{tj)l/ " ClP])klm " Opklmn 



3/ n'-'i,Y„^V'.«¥f.Y«^) (3.11) 



V Vh dpj)ki Vh dpk{^j) J 

up to surface integrals which we omit. 

We have checked that, as they must, the two constraint equations (|3.7p with p.Sp reproduce the (00) and the (Oi) 
components of the vacuum Euler-Lagrange equation of motion (|1.2p and that (j3.9|) and (|3.10p are nothing but their 
{ij) components modulo constraints (the calculation is fairly involved and the details will be presented elsewhere ,13J). 

To compute the Poisson brackets of the secondary constraints p.Sp , it is useful to define smeared quantities by 

H[u] = [ d"-^xCi^, M[C]^ [ d^-^xC^C, (3.12) 

where v{x^) and ^'^{x^) are test functions. This allows us to integrate by parts and to ignore boundary terms, 
as everywhere else in this paper. A calculation shows that the Poisson brackets of the smeared Hamiltonian and 
momentum constraints, modulo p.8p . read the same as in General Relativity: 

{M[^ilM[Q} = M[CiD,Q - e^D.a] , (3-13) 
{M[e],H[iy]}^H[£^iy]. 

Here it may be worth making a few comments on the relation between the constraint p.Sp and the Hamiltonian and 
momentum constraints. 

If p.Sp is invertible with respect to O^j , we may eliminate Vlij from the action completely, leaving no further 
constraints. Then the remaining constraints are the Hamiltonian and momentum constraints, which are of "first- 
class" [l^HJl representing the diffeomorphism invariance of the /(Riemann) action. 

On the other hand, if p.8p is not completely invertible, some of the components give rise to non-trivial extra 
"primary" constraints on the dynamical variable 5'*-' . Then time derivatives of these extra constraints may give rise 
to "secondary" constraints 23]. After all the extra constraints, irrespective of primary or secondary, are spelled out, 
one can classify them into first-class and second-class. In most cases, thesj3 constraints will be of second-class, which 
may be inserted in the action to reduce the dynamical degrees of freedom [16] . Then, for consistency, the Hamiltonian 
and momentum constraints expressed in terms of the reduced phase space variables should also satisfy Eq. p.l3p . As 
examples, this will be confirmed below in the case of Einstein gravity as well as of /{TV) gravity. 

In some exceptional cases, these extra constraints may happen to be of first-class, that is, the action may acquire 
a larger gauge invariance: It is known that the "Weyl^" action in Z? = 4 has a conformal invariance and in that case 
the constraint algebra of the Poisson brackets is extended to incorporate the generator of conformal transformations 

Now that we have completed the presentation of the general formalism, let us turn to some "practical" applications. 



IV. PHASE SPACE REDUCTION 



In this Section, we show how the constraint p.8p is utilised to find "reduced" Hamiltonians for various sub-classes 
of /(Riemann) gravity. We follow the procedure advocated in ^1^: Second-class constraints arising from (|3.8p are 
inserted into the first-order action p.2p to eliminate as many components of flij and of dynamical variables as 
possible. If there remain any constraints on ^f*^ , irrespective of whether they give rise to further constraints on the 
other dynamical variables or not, it becomes necessary to redefine canonical momenta conjugate to the reduced sets 
of hij and '^'^^ . As we shall see below, they are relatively easily read off from the first-order action. Consequently the 
Hamiltonian in terms of the reduced set of variables is obtained. 
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A. / = 7^^.p<. w'""' 

This is a simple "generic" example where all the components of Vlij can be extracted from p.Sp , so that there are 
no extra constraint on the dynamical variables. 
The first step is to decompose f{Q^upa) as 

f{Qt.vpa) = Q^upa (f'""' = 4 • + p^u p''^' - 4 p.,fe p''"" , (4.1) 

where pijia and pm^ are given in (j3.6p . Thus the constraint (|3.8p reads 

^''=-—- (4-2) 

This constraint is inserted into the first-order action (|3.2I) to eliminate . The Hamiltonian density is therefore 
given by l|3.4p and (13. 5p . where 

* • + ^ = - " + \ {P^,M P^^'' - 4p,,fe p^^^) . (4.3) 

It depends on the usual set of variables of General Relativity, {hijjp^^ , N, P^} , plus D {D — l)/2 extra degrees of 

freedom, { ^-^^ , } . 

Let us mention that the form of the canonical equations of motion (|3.9p and (|3.10p is unchanged even if (|3.8p is 
inserted into the Hamiltonian before taking the variations. This is because the Hamiltonian depends on fiy only 
through the combination 

^■n+^f{g^,p^), (4.4) 

so that the additional terms appearing in the equations of motion due to the constraint (j3.8p are all proportional to 
the derivative of (|4.4p with respect to fiy , which vanish by virtue of p.Sp itself. 

It is also worth mentioning here a decomposition of \E''-' in the "generic" case. The D {D — l)/2 components of 'J*-' 
can be decomposed into two irreducible parts: 

$ = -^^, v^^^ = T**^ (4.5) 

where ^ = hij ^'^^ and the symbol T denotes the traceless part. Noting that there is an arbitrariness in choosing 
canonical momenta conjugate to the variables {hij, V'*"'} i one finds it most convenient to introduce 

= P'' - (n + ^ ) , TT,^ = tH„ (4.6) 

resulting in 

C* =p- h + n^ + TT -ijj ^n* , (4.7) 

where Ti.* is to be given in terms of the new variables. Now it is manifest that the "generic" theory contains a scalar 
degree of freedom {$, H} and traceless tensor degrees of freedom {V'*"', T^ij} having {D + 1) {D — 2)/2 components on 
top of the canonical metric degrees of freedom {hij,p'^} . 

B. f = n 

We show here how the ADM Hamiltonian for General Relativity follows from our general formalism. 
The decomposition of f{Qpvpa) is 

f = Q=crg''''Q,,,p„^-2^ + p, (4.8) 
where 57 = K^^ ftij and p = K^^ h^' pijki so that the primary constraint (|3.8p reads 

= h'^ . (4.9) 
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Hence it "freezes" all the extra degrees of freedom of ^'^-^ . Though (|4.9p does not give Jly , the terms containing fiy 
in the Hamiltonian will all drop out due to the linear dependence of / on 0^ . 

The Hamilton equations l|3.9p for the velocities give a secondary constraint. It tells us that, since ^t''-' = h''^ , the 
momentum conjugate to 4'*-' is also frozen: 

n,j^2p,j-^h,j, (4.10) 

where p = hij p^^ . 

We now follow the procedure advocated in [l^ and insert (|4.9p into the first-order action (13. 2|) , which becomes 
(ignoring the divergence) 

C* ^p-h-n* , (4.11) 

with 

f^=-f^ + ^h'^. (4.12) 

This p''^ plays the role of the new momentum conjugate to hij . 

We now gather the results, to wit, all the dynamical variables are expressed in terms of {hij,p^^} as 

Mj^^=h^=^ n,, ^-2p,, + -^/i,, , ^-f^ + ^^h^^ , (4.13) 



where p = hij j?^ , and we plug them into the Hamiltonian (13. 4p and ()3.5p to obtain 

n*=NC + l3'Ci, (4.14) 

with 



(4.15) 



This is nothing but the ADM Hamiltonian for General Relati vity in D dimensions. As for the equations of motion 
(13. 9p and (|3.10p . they reduce to the ADM equations, see e.g. [l^l. Moreover, the constraints (|4.15p are first-class as 
is well known. 



c. / = /(7^) 

We show here how our general formalism yields the "Jordan- frame" Hamiltonian of f{TV) gravity. 
As above we have g = —2 Vl + p with Vl = /i*-' ilij and p = pijki , that is, using p.6p : 

H^ - H • H 

Q = -2n+ + (4.16) 

h 

so that Eq. (|3.8p reads 

^^^^f'{e)w^ ^ = ^vi/'i^o, (4.17) 

where the symbol T denotes the traceless part. This primary constraint tells us, first, that the trace 17 is known 
in terms of $ and other variables, and that the traceless part of ^''^ is constrained to be zero so that ^''^ reduces 
to one scalar degree of freedom, $ . The traceless part of $7^ automatically disappears from the Hamiltonian since 
^ff • i7 = $ $7 and f^g) depends only on the trace from the beginning. 

The traceless part of Hy is also constrained through the Hamilton equations for the velocities p.9p . Indeed, the 
traceless part of the velocities now satisfies t(^'"') = i{h^-') giving a secondary constraint 



Tlljj = — TPrj 



(4.18) 
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In order now to find a new momentum p^^ conjugate to hij , we again follow ^Td\ and insert (|4.17p and (|4.18p into 
the first-order action (j3.2p . which now reads (ignoring the divergence) 



C* =p - h + Il<i>-H* , 



(4.19) 



with 



D- 1 



(2p-<i>n), 



(4.20) 



where p = hij p^^ . Gathering the results: 



-2 p.,, 



1 



(2p + $n) 



pv ^ _pV 



where p = hij p^^ , and using (|4.16p and (I4.17[) , we have 



h^^ 
D-1 



(2p + $n) 



^'i>[U + 2 (i? - 1) /i $ h 



(4.21) 



(4.22) 



where it is understood that g is known in terms of (f> via f'{g) = ^. Therefore the Hamiltonian (13. 4|) becomes 
n* =NC + P'C, with 



Vh 



TP ■ TP 



D-2 



pU 



D- 1 



(4.23) 



The constraints (|4.23p give the standard Jordan-frame Hamiltonian for f(jV) gravity, see Q. Moreover, it can be 
checked that they are first-class. 

One can also make a canonical transformation of the Hamiltonian into that of Einstein gravity with a minimally 
coupled scalar field. These metrics are related by a ZJ-dimensional conformal transformation, see e.g. [2^] and 



Df — r ptJ,vp<T 

We consider here the case when / = Cf^i^p^ C'"''"^ where C^^pcr is the I?-dimensional Weyl tensor. In terms of the 
Riemann and Ricci tensors, it is expressed as 

/(7^^.p.) = 7^p.p. n^"""- - n,, n^" + ■ (4-24) 

This theory yields conformally invariant equations of motion in Z3 = 4 . 
To decompose f^Q^vpa) , we need 

e^'' = r^^ + r2 • r2 - 2 p • r2 p • p - 2 (4.25) 
^^2 = (-2^! + p)^ 

where = Qij , py = /i''' pikji , p = h'-'' p^jki and pi = h^'' pjik with 



H.fcH,,-H.,H,-fe HH,,--(H.H),, - n ■ H 
Pijfci — h Kijkl , Pij ~ H -Kij , p — h it , 

For comparison with the Riemann squared case, see (|4.ip and p.6p . The constraint equation (j3.8p now reads 



(4.26) 



- - [{D - 3) T^^'^^' + TP''] , (4.27) 
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where again T denotes the traceless part. Thus only the traceless part of i^ij is determined. As a consequence, there 
appears an extra primary constraint h ■ = , that is, the trace part of 5'*-' disappears from the action. This imphes 
the number of degrees of freedom in ^P*-' will be reduced by one. 

Now we have a relation between the velocities h-'^ + h- '^ = Oasa consequence of the above primary constraint 
so that the Hamilton equations for the velocities (|3.9p yield a secondary constraint. 



(4.28) 



D 

T*-Tn = o, 



where p = hij p''^ . In contradistinction with the previous example of f{TZ) gravity, this equation does not constrain 
n but p. For the moment, we decide not to insert this constraint into the action. The reason will be clarified below. 



In order now to find the new momenta p^^ and tt^ conjugate to hij and ^p"^- 
[l6t and write the first-order action (|3.2p as (ignoring the divergence) 

C*=p-h + T:-ip-n\ 



with 



Let us gather the results: 



p^3 



n,: 



n 



D - 1 



n 



D- 1 



p^3 



tVE'*-' , respectively, we again follow 

(4.29) 

(4.30) 



n 



D - 1 



where we see that 11 is not determined. Plugging these into the Hamiltonian p.4p one finds 

2 ivn 



7VC + /3*Cj + 



with 



1 



C = — = ( — tp ■ n ■ TT — 2p ■ n] + Vh 



D-1 Vh 
D-2 

D-3 ■■ " ■■ J ■ • [8(1) -3) 

TPijkl TP'^'"^ - 4:TPijk TP'^'') , 

TT^kD.i^^'' + 2VhDk 



Cw , 



(4.31) 



(4.32) 



D-2, 



Vh 



(4.33) 



Cw = 



where p = hij p*-' . We note that the Hamiltonian is linear in H . In other words, H acts as a Lagrange multiplier 
whose equation of motion, Cw = , coincides with the secondary constraint (|4.28p . 

Boulware @ showed that, in _D = 4, the Poisson bracket algebra of the constraints closes, that is, Cw as well as C 
and Ci are first-class constraints. Boulware also showed that Cw is the generator of conformal transformations under 
which the theory is invariant. Thus, to summarise, the Lagrangian for the Weyl squared theory in D = 4 keeps the 
original form (|4.29p . 



C* =p-h + T: -iP-n* , 



(4.34) 



where TL* is given by (|4.32p . but with all the constraints C, Ci and Cw being first-class. 

In passing, we note that it may be more transparent to separate out the determinant of hij as an independent 
canonical variable. Namely, in addition to p= hij p^^ , introducing 



r; = ln/ii/(^-i) 



h 



e" TP'^ 



(4.35) 



we express the Lagrangian as 



C* = p- h + pr] + TT-xp — Ti* , 



(4.36) 
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where 7i* is given by 

n* =NC + I3' C; + WCw (4.37) 

Here we have replaced nbyW^=2(_D— 1)^^A^ VTi 11 , and C , Ci and Cw are to be expressed in terms of the new 
canonical variables {/ly , p*-'}, {f7,p} and TTij} . 

In > 4, however, the Poisson brackets oiCw with C and no longer close, giving rise to additional secondary 
constraints. This is because, although the Z?-dimensional Weyl tensor, C^^pc, , is conformally invariant in any dimen- 
sions, the action is not conformally invariant any longer in dimensions other than D = A. In this case, there may 
be more secondary constraints from the time derivatives of these secondary constraints. They will be all second-class 
constraints. Inserting them into the action will reduce the phase space considerably. Because of rather involved 
calculations, we have not checked how many second-class constraints would appear in the end. We plan to come back 
to this issue in a future publication [isj . 



V. CONCLUSION 



We have presented in this paper a canonical formulation of /(Riemann) theories in a form as compact as possible. 
It includes in a unifying manner well-known subcases and should prove useful to analyse the properties of various 
theories of "extended gravity," such as the global charges associated with the solutions, their (in)stability and the 
(non-)positivity of energy or the jmiction conditions. 
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APPENDIX A: THE OSTROGRADSKY HAMILTONIAN OF /(Riemann) GRAVITY 

For completeness we give here the canonical transformation which relates the variables used in this paper to the 
"Ostrogradsky" one, that is, the extrinsic curvature of the ADM foliation, see Q for an Ostrogradsky treatment of 
quadratic theories. 

The Ostrogradsky action written in terms of <r^,ypCT = 'T^^lupa\Kij~*Qij , where Qij is independent of other fields, is 



S= / d^a;V=^ 

IM 



(Al) 



The momenta canonically conjugate to hij and Qij , respectively, are 



fc^.^i^ = V7Iu^ P^=^1^ = -^1L^ (A2) 
where T,ij = n'^ ^i^u ■ We have the following relations with the variables of this paper: 

ny = -V7^g.y, (A3) 



They give 



p-h + n--^ ^k-h + P -Q- ^^{Q ■ P) (A4) 
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with 

fcy ^ pv + /jy" . (A5) 
Plugging these relations into our Hamiltonian (|3.4p , we obtain the Ostrogradsky Hamiltonian H* ^ N C + P'^ Ci with 



^ (A6) 

where '^fj.upa- is expressed in terms of hij , Qij and , and where S^j is constrained to satisfy 

=-2^. (A7) 
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